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1 Introduction 

Propagation of sound through the turbulent atmosphere is a statistical problem. The randomness of the 
refractive index field causes sound pressure fluctuations. Although no general theory to predict sound 
pressure statistics from given refractive index statistics exists, there are several approximate solutions to 
the problem. The most common approximation is the parabolic equation method. Results obtained by 
this method are restricted to small refractive index fluctuations and to small wave lengths. While the first 
condition is generally met in the atmosphere, it is desirable to overcome the second. This paper presents 
a generalization of the parabolic equation method with respect to the small wave length restriction. 

2 Parabolic Equation Method 

For the small wave length limit 1 the Helmholtz equation can be converted into a parabolic form (main 
propagation direction e z ) f\/: 

( 2 ^' + S? + & + *'■«) 0M = O M 

k = wave number; p = refractive index deviation; ip = complex sound pressure 

The refractive index deviation p is considered as a random function. Therefore equation (l) is a stochastic 
differential equation and the sound pressure tp becomes a random function, too. Stochastically this equation 
is nonlinear, e.g. it contains a product of two random variables. For this reason the stochastic parabolic 
equation cannot be solved exactly. Further approximations are necessary. Several mathematical tools 
were applied to remove the stochastical nonlinearity, i.e. path integrals /2/, functional derivatives /3/, 
perturbation expansion methods /!/. Despite approximations used in the calculations looking different, 
the results are all the same. 

The physical meaning of all these approximations becomes evident in the local method of small perturbations 
/ 1/ . By this method the scattering volume is divided into slabs perpendicular to the main propagation 
direction. Each slab is chosen as thin as required by the validity limit of the first order perturbation 
expansion term (single scattering approximation, Born approximation). This distance clearly depends on 
the strength of the refractive index fluctuations. If these fluctuations are sufficiently small, the slabs are 
much thicker than one correlation length of the random medium. Therefore the slabs can be regarded as 
uncorrelated. Based on both assumptions - small refractive index fluctuations and a small wave length 
compared to the correlation length - the statistical independence of subsequent slabs can be proofed 
mathematically. Wave propagation through random media is described here as a Markov process. 

*The wave length must be small compared to the size of a typical inhomogeneity of the medium. In statistical terms this 
size is expressed by the correlation length of the refractive index autocorrelation- function. 
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As a consequence of the Markov property slabs of finite thickness are no longer necessary. This results in 
a differential equation for the mean sound pressure which is linear in the stochastical sense as well 2 /!/: 

< t/>(r) >= 0 (2) 

( 3 ) 


< V'(r) >= ^o(r) exp{-az} (4) 

tpo — free propagated incident wave 

The mean complex sound pressure decreases exponentially while the wave propagates through the random 
medium. This is an effect of decorrelation of the sound wave due to phase fluctuations. Different members 
of the statistical assembly interfere destructively because of their different phases. 

The validity of result (4) - and of any other result obtained by the parabolic equation method - is restricted 
first by the validity of the parabolic wave equation (1) and second by the validity of the Markov assumption. 
Necessary conditions are the smallness of refractive index fluctuations and the smallness of the wave length 
(compared with the correlation length). In the next section the first condition is also assumed to be 
true. The small wave length assumption, however, will be replaced by the weaker condition of negligible 
backscattering. This will lead to a generalized Markov process and, consequently, to generalized results 
with respect to the wave length - correlation length ratio. 

3 Generalized Local Method of Small Perturbations 

While generalizing the parabolic equation method the main idea of the local method of small perturbations 
will be retained. The refractive index fluctuations are assumed to be small enough to justify the application 
of a single scattering approximation within a distance A z in the scattering volume which is large compared 
to the correlation length. Again the scattering volume is divided into slabs of this size. Therefore subsequent 
slabs are uncorrelated as well. But contrary to the parabolic equation method the contribution of one slab 
will not be calculated from the parabolic equation but from the Helmholtz equation: 

(a + k 2 ( 1 + M (f))) tp(r) = 0 (5) 

A = Laplace operator 

Neglecting backscattering 3 yields a difference equation for the mean sound pressure: 

< ip(p n ,nAz) >= (g 0 + £) < (n - l)Az) > (6) 

2 It is possible to derive similar equations for the higher statistical moments of V* by the parabolic equation method, too. 
Only the first moment equation and its solution are presented here in order to compare them with the generalized results 
derived in section 3 of this paper. The decorrelation coefficient a is calculated for an exponential autocorrelation function for 
comparison, too. 

3 Neglecting of backscattering is not an assumption of the parabolic equation method but a consequence of the small wave 
length limit. The parabolic equation implies that there are small scattering angles only and no backscattering at all. 


/ d d 2 d 2 \ 

( al + + v + j 


a = 


<n 2 > k 2 i 


l = correlation length 
Equation (2) can be solved easily: 
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p — (x, y) ; n = number of slab 

A A 

Go is the integral operator of the homogeneous Helmholtz equation (p = 0, free propagation) and S is 
an integral operator for the scattering within one slab. Since double scattering is the lowest order non- 

A 

vanishing term, the kernel of S contains the autocorrelation function of the refractive index field: 

S {r nt r) < tp(f) >= -k 4 Jd*? Jd?r G(f n , f')G(f' ,r) < > < H*) > ( 7 ) 

G = Greens function of the Helmholtz equation 

The solving of equation (6) is somewhat different from that of the related parabolic equation problem. 
Here no wave length approximation helps to calculate the integrals. But for the case of a homogeneous 
refractive index autocorrelation function, equation (7) becomes a convolution product. It can be Fourier- 

transformed with respect to the variable p, and this operation turns the operators S and Go into simple 

ni rs/ 

functions S and Go. In the Fourier representation equation (6) reads: 

<V> (k, nAz) >— + s^j <ip ( K , (n - l)Az) > (8) 

k = 2-dimensional spatial frequency 

By Fourier-transformation the ^-integrations are already performed and only the z-integrations are left. 
They can be performed, too, if the z-dependence of the medium autocorrelation function is known. For 
the sake of simplicity this dependence is assumed to be exponential 4 . After z-integration the scattering 
contribution of one slab is seen to be proportional to the slabs thickness Az. Therefore equation (8) can 
be written as (5 = s Az): 


<V» (/c,nAz) >= + 8 A <«/> (/c, (n — l)Az) > 

The effect of all slabs is obtained by iteration: 

<^> (k, nAz) >= ^G 0 + s A z^ (k,0) 

Regarding a sufficiently large number of slabs yields (z = nAz): 

»v /v 

<t/> ( k,z ) >=V> 0 («.*) exp{a (/c)z} 

For the special case of an isotropic exponential autocorrelation function s becomes: 


( 9 ) 


( 10 ) 


( 11 ) 


< ft 2 > k 4 
4 oaj(aj - a) 


( 12 ) 


o = Vk 2 — /c J , aj = -^(k + i/l) 2 - k 2 (13) 

To compare (11) and (12) with the corresponding parabolic equation method results (4) and (3), the incident 
wave V»o is assumed to be a plane one. Then the Fourier-transformation of (11) results in: 

4 For more complicated functions the s-integration can be performed numerically. 
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(14) 


< i K*) >= V>o(0 ex p{ * (°) z ) 


5(0) = 


< m 2 > * 2 / (i« - J fe 2 / 2 ) 

4 (1 + fc 2 / 2 ) 


(15) 


The real part of 5 (0) describes the decorrelation of the sound wave. It is a more general expression than 
(3) - only in the small wave length limit they are equal. Because of s(0) being a complex number, a second 
effect is predicted by this method, which cannot be seen in the parabolic results. The imaginary part 
is a stochastic correction to the wave number k due to an increase of the mean propagation distance in 
the random medium. Only in the small wave length limit the scattering angles are small and the mean 
propagation distance corresponds to the z-extension of the scattering volume. 


4 Conclusions 

A generalized form of the local method of small perturbations has been presented in this paper. Working 
directly from the Helmholtz equation instead of the parabolic equation the small angle scattering method 
was replaced by a forward scattering method . By this method only one result wets derived here: The first 
statistical moment of an incident plane wave scattered by a very weakly statistical homogeneous random 
medium with an exponential autocorrelation function. This result shows corrections to the corresponding 
parabolic equation method result. 

It is possible to apply the method to more complicated problems, i.e. a difference equation for the second 
statistical moment can be derived by the same idea. 

If the medium fluctuations become stronger, the thickness of one slab decreases. The slabs might be thicker 
than the correlation length, but not as much as assumed before. Then correlations between two slabs have 
to be taken into account. This leads to difference equations which connect statistical moments not only 
from one slab to the following, but to the next following, too. 

The most valuable advantage of this method might be its suitability for numerical calculations. For any gi- 
ven medium autocorrelation function the scattering function s can be obtained by Fourier-transformation. 
The incident wave is also Fourier-transformed. Then the algorithm given by equation (9) is applied iterati- 
vely until the desired propagation distance is covered. The final result is obtained by Fourier-transformation 
again. 

The author wishes to thank Prof. K. Haubold, Prof. V. Mellert, Dr. M. Schultz-von Glahn and A. Sill for 
the valuable discussions during the whole work. 
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